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Abstract 

We show that the analytic single-particle density of states and the optical 
conductivity for the half-filled Hubbard model on the Bethe lattice in infinite 
dimensions describe quantitatively the behavior of the gap and the kinetic 
energy ratio of the correlated insulator V2O3. The form of the optical con- 
ductivity shows a; 3 / 2 rising and is quite similar to the experimental data, and 

the density of states shows u 1 / 2 behavior near the band edges. 
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It has been an interesting question that how well does the Hubbard model describe a real 
strongly correlated system. To answer this question, appropriate experiments are required. 
Recently, Thomas, et. al. |1[ have measured optical conductivity for the insulating state of 
V2O3 whose magnetic phase has been known as antiferromagnetic 0. They observed some 
interesting features for the properties of the insulating phase of V2O3 from the measurement 
of optical conductivity. They found that the insulator gap is not the Slater gap come from 
antiferromagnetic ordering || but the correlation gap resulted from strong on-site repulsion 
d]. In other words, the optical conductivity rises as {u — 2A) 3//2 , where 2 A is the gap 
width, instead of to 1 / 2 expected from a Slater antiferromagnet. This uj 3 ! 2 rising is expected 
from the to 1 ' 2 behavior of the single-particle density of states (DOS) near the band edges. 
They compared their experimental data with various theories |^-||. The best fitting one 
was the paramagnetic solution of the dynamical mean-field theory 0|| which is valid in 
infinite dimensions. They explained that the reason why paramagnetic solution describes 
the insulating state of V 2 3 quite well was caused by the spin frustration. The spin structure 
of V2O3 is substantially frustrated, i.e., there are both ferromagnetic and antiferromagnetic 
nearest neighbors. In addition, it is known that a theory in infinite dimensions is well- 
applicable to a bulk system. |J They also argued that their results are representative of the 
zero-temperature limit. 

In this Letter, we report that the behaviors of uj 3 / 2 in optical conductivity and uj 1 / 2 
in single-particle DOS can be seen in the Hubbard model on the Bethe lattice in infinite 
dimensions and our results in infinite dimensions fit experiments over all the insulator regime 
of V2O3, while the dynamical mean- field theory shown in Ref. 1 does not cover full insulator 
regime, especially near the metal-insulator transition. We show explicitly that the form of 
the optical conductivity is quite similar to that of experiment. In addition, the area under 
the curve of optical conductivity which gives the average kinetic energy has been compared 
with experiment by using the ratios of the averaged kinetic energy to its noninteracting value. 
The agreement between our theoretical values and experimental data is quantitatively good. 

We now show our theoretical works and comparison with experiment. Before we go 
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further, we first consider spin frustration in the Bethe lattice since V2O3 is considered as a 
spin frustrated system . We include next nearest neighbor hopping to take the frustration 
into account, since the Bethe lattice is a bipartite lattice which does not have frustration in 
lattice arrangement. For the Bethe lattice, however, the number of next nearest neighbor 
is q(q — 1), where q is the coordination number. Therefore, t 2 representing next nearest 
neighbor hopping integral must be scaled as ti = t*/q to make the kinetic energy finite in 
infinite dimensions. The t 2 , however, is 1/y/q times less than t\ which represents nearest 
neighbor hopping and scales as t\ = U/y/2q flOfl . Thus the effect of next nearest neighbor 
hopping can be neglected in calculating the DOS on the Bethe lattice in infinite dimensions. 



Therefore, the result obtained in the previous work |TT[ can be used without any change for 
the present analysis. 

We briefly introduce the dynamical Lanczos method P, p!2Hl3| used in this work. The 
single-particle DOS can be obtained by calculating the one-particle Green's function of 
the fermion operator at the same site, i.e., (^o|{ c J CT (£)> Cjo-ll^o) where cj CT and Cj a are the 
fermion creation and annihilation operators with spin o at site j, the curly brackets mean 
anticommutator, and |^ ) denotes ground state. The single-particle DOS p a {uj) is given by 



p*W = £ ImG §"V + «0> (i) 

where 

= ^(^o|{c JCT , (u + L + ie)-%}\* ) 

= — ^jj{z)\z=-iL>+ei (2) 

where L is the Liouville operator. The superscript (+) denotes the usual notation of the 
retarded Green's function [|J. The on-site Green's function ^jj(z) can be represented by an 



infinite continued fraction 12.13 
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z - a + 



(3) 



z— ai + - 



where a v = (iLf v , /„)/(/„, /„), A^ = (U, fu)/(fu-i, fu-i)- The inner product is denned 
by (A,B) = (^q\{A, -B^jl^o); where A and I? are operators of the Liouville space, 

fit 

is the adjoint of B. These a v and A„ are obtained using a recurrence relation /„ +1 = 
iLf v - a v f v + A„/„_i. [0 

The Hubbard model which we study here is written as 

U 



H 
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By choosing / = Cj a , we have obtain the on-site Green's function Ejj(z) for the half-filled 
Hubbard model on the Bethe lattice in infinite dimensions whose ground state is assumed 



as paramagnetic, [fUJl 
where A 1 = f + § , A 2 = (U 2 + t 2 )/(4A 
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If we set the chemical potential at /i = Eq. (5) gives the single-particle DOS for the 
insulating phase (a > b) as follows: 



-ReZjjiz) 
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-iuj+0+ 
AiA 2 
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A 2 



(6) 



f (6-a) + A 1 + (f -IV 2 

where H 7 = {cu 2 - (-y/a - Vb) 2 }{(\/a + v^) 2 - ^ 2 }- We take (-) sign for u > and (+) 
for uj < to satisfy the boundary condition Hj 3 -(t = 0) = 1 given in Eq. (2). The lower and 
upper Hubbard bands exist — (\fa + \/b) < ^ < — \fa. + Vb and y/a — Vb < u < \fa + Vb, 
respectively. One can observe uj 1 ! 2 behavior near band edges from Eq. (6). 

One can observe that this single-particle DOS has band width 2D = 2\/b = 2t* and band 
gap 2 A = 2(\/a — Vb) = U — 2D. Thus we get the following relation: 
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We draw Eq. (7) with experimental data in Fig. 1. A remarkable agreement is seen over all 
insulator regime. 

We now obtain the optical conductivity cr(cj) using a formula valid in infinite dimensions 
0, 



a(co) = (T / dJ I dep® (e)p(e, u')p{e, J + uj) f ^ + ^ , (8) 



UJ 



where f(uj) is the Fermi distribution function and o"o = nt *2%v — where a, N, V are lattice 



constant, number of lattice sites, volume, respectively. In Eq. (8), p^°\e) = ^J2t 2 — e 2 
which is the single-particle DOS for U — 0, and p(e,uj) = — -haG(e, uj) = —^Im{u + irj — 
e — S(cj)] -1 . Use of the momentum-independence of the self-energy in infinite dimensions 
has been made. This property make it possible to express the self-energy in terms of the 
on-site Green's function G{uS) = —iEjj(—iuj). 

If we set ( = uj — £(o>), the one-particle Green's function is mapped into the frequency 
renormalized noninteracting one which describes the noninteracting system under effective 
field, i.e., G (0) (C,e) = [C - e] _1 - Since we obtain a v = iU/2, A v = t 2 j2 for all v for the 
noninteracting Hubbard model on the Bethe lattice in the paramagnetic state, Eq. (3) gives 



w (o)/- : \lz 2 + 2tl 



'33 ^ t 2 4 t 2 



= iG (0 \() = iG{uj). (9) 

-iC+0+ 



Solving Eq. (9) for ( gives the self-energy as 



E(w )= u -| G(u) - g ^. (10) 

This relation has been obtained by Georges and Krauth |7| in terms of the effective action 
theory. 

Using Eqs. (5) and (10), we get p(e, u) and finally the optical conductivity a(u) from Eq. 
(8). Fig. 2 shows theoretical a(uj/t*) in units of <To for various U/t*. We use approximation 
Ai pa a and A 2 pa b in drawing a(u/t*). We choose U/D = 2.1 and U/D = 4 to compare 



with experiment. To make the comparison appropriate, we need to adjust ou/t* for each 
sample, since each sample has different band width D which is equal to t* in our theory. 
Therefore, for the horizontal scale, we set one unit of u/t* to 0.31eV for U/D = 4 and 1.21eV 
for U/D = 2.1. The former gives 2A = 0.62eV and the latter 0.12eV. Both gap widths are 
within experimental error bounds. For the vertical scale, however, we multiply each ct(uj) 
by the corresponding D^ (D m = 0.47eV and 0.31eV for U/D = 2.1 and 4) in Ref. 1 to take 
o"o which is proportional to t\ into account. Since the vertical scale itself is arbitrary, only 
relative height is meaningful. The optical conductivities represented by experimental scale 
are shown in Fig. 3 with lu 3 / 2 rising expressed by the dashed line. 

Even though the theoretical value of the optical conductivity at a particular frequency 
is not quite close to that of experiment, overall structure is quite similar each other. An 
additional test for our theory related to optical conductivity can be performed by comparing 
the ratio of the average kinetic energy (T) to its noninteracting counterpart (T)o with 
experimental data. Since the conductivity sum rule 



gives rise to the average kinetic energy, one can immediately obtain the average kinetic 
energy by performing the integration. The explicit form of the constant £ has been given 
by £ = jj^z [[y|. However, we do not need the explicit expression, because it is cancelled in 
getting the kinetic energy ratio. 



in Eq. (11) for the optical conductivities shown in Fig. 2, we obtain the ratios as follows: 
(f )/(f ) = 0.349, 0.205, 0.147, and 0.116 for U/D= 2.1, 3, 4, and 5, respectively. We put 
these values with experimental data and other theoretical work in Fig. 4. Our theoretical 
values are quite close to the result of dynamical mean field theory and its extension to the 



In conclusion, we argue that the paramagnetic solution obtained by the dynamical Lane- 




(11) 




. Using this and performing the integration 



metastable branch where both insulating and metallic solutions coexist 
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zos method for the Hubbard model on a Bethe lattice in infinite dimensions may describe 
the insulating phase of the three-dimensional strongly correlated system like V 2 3 quantita- 
tively well. The spin frustration effect which makes the system paramagnetic-like has been 
treated by considering next nearest neighbor hoppings in the Bethe lattice. The behavior of 
the insulating gap and the average kinetic energy ratio according to the change of U/D is 
quite well agreed with experiment quantitatively, while the optical conductivities showing 
(uo — 2A) 3 / 2 rising are agreed qualitatively. Finally we expect that the metallic regime of 
V2O3 can also be explained quantitatively in terms of our result for metallic phase. This 
will be a forthcoming work. 
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Figure Captions 



Fig. 1 : Comparison of Eq. (7) with experimental data (solid circles for optical measure- 
ments, open circles and diamonds for dc measurements) of Ref. 1. The dashed line is 
the present theory and the solid line is the stable solution of the dynamical mean field 
theory. 

Fig. 2 : Theoretical optical conductivities for U/D — 2.1, 3, 4, and 5 obtained by using Eq. 
(8). We express a(cu/U) in units of a . 

Fig. 3 : Comparison of optical conductivities for U/D = 2.1 and 4 with experimental data 
(solid points for U/D ~ 2.1 and open circles for U/D ~ 4) shown in Ref. 1. The 
solid lines are theoretical values and the dashed and the dotted lines denote u; 3//2 rising 
in theory and experiment, respectively. The horizontal scale for theoretical curves is 
reexpressed by the energy scale used in experiment. Arbitrary units are used for the 
vertical scale for theoretical values. Details are explained in the text. 

Fig. 4 : Comparison of average kinetic energy ratios (T)/(T) for U/D = 2.1,3,4, and 5 
with experiment (solid circles) and theory. Present theoretical values are expressed 
by solid triangles. The solid, dashed, and dash-dot line denote dynamical mean field 
theory given in Ref. 1. The region covered by dashed line is the metastable regime 
suggested in Refs. 7 and 8. 
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